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Abstract 

We explore the feasibility of estimating primary cosmic ray composition at ultra high energies from the sum of 
muon, electron and photon densities and the depth of maximum of extensive air showers detected by the Auger 
Observatory. From the information of the Fluorescence (X max ) and water Cerenkov detectors (p(1000)) 
we infere the most probable type of primary which originated the shower. The method is tested simulating 
Extensive Air Showers (EAS) at energies up to 10 20 eV using a version of MOCCA/Sibyll adapted for Dec 
Alpha Servers. We also discuss some results on the mass composition using the depth of the shower maximum 
j and the lateral distribution of showers in the energy range 10 18 -10 20 eV. 
OV 

^ ■ 1 Introduction: 

In recent years ultra-high-energy cosmic rays have received renewed attention from the whole physics com- 
munity. Most of the interest for this subject comes from the observation of eight extremely energetic particles 
(Auger Collaboration, 1997), inspite of the so-called GZK cutoff. In 1966 Greisen and, independently, Zat- 
sepin and Kuzmin predicted the suppression of cosmic rays of energies as high as W 20 eV or above, an effect 
£Nj \ known as the GZK cutoff. As they pointed out, this suppression is simply due to the fact that protons with 
energy above about 4.1 x 10 19 eV should interact with the cosmic microwave background radiation, implying 
that no such high-energetic particles should be observed. 
, Since there had been observations of those air showers that seem to have been induced by cosmic rays of 

that high energies, the intriguing challenge is that possible sources of these particles should not be located far 
^ ■ away, otherwise they would have been suppressed by the GZK effect. The details about those events can be 
O . found in the Design Report of the Auger Project (1997). 

In order to observe air showers induced by UHE cosmic rays, the Auger Collaboration has chosen to com- 
bine two different detection techniques already explored by two long-lasting collaborations. These techniques 
involve the observation of the Cerenkov light emitted by ultra-relativistic electrons and muons in water tanks, 
as used in the experiment of Haverah Park (Lawrence, Reid & Watson 1991), and the observation of the fluo- 
rescent light with mirrors and photomultipliers, as done in the Fly's Eye experiment (Baltrusaitis et al 1985). 
_C ■ This fluorescent light is emitted by N% molecules in the atmosphere when excited by air shower particles, 
J> ! mainly electrons and positrons. 

I Measurements of arrival times give an estimate the shower direction, also assumed as that of the primary. 

From this information and the signals in the water tanks and/or fluorescent detectors, the fraction of the primary 
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energy gone into electromagnetic particles and the depth of the shower maximum can be inferred. Those 
combined observations will give an estimate of the primary mass composition and also a hint on the possible 
origin of those extremely energetic particles. 

Unfortunately the measurements of both Haverah Park and Fly's Eye experiments do not allow to estimate 
the mass of the primary at an event-to-event basis, but only a mean composition for a given energy range. 
Sokolsky et al (1992) presented a detailed report on this subject, to serve as a reference guide about the 
detection techniques, their experimental results and methods to characterize UHE primary cosmic rays. From 
this report, the difficulties inherent to this intend are clearly seen. 

Basically the quantities which have been used to estimate primary mass composition are: 
1) The ratio between muon and electromagnetic densities; 



I) Measurements of the depth of shower maximum (X max ) and its distribution; 

3) The elongation rate (ER), or the mean rate of increase of X max with primary energy (Eq). 

We have tried to explore some correlations between those quantities for simulated air showers, in order to 
find out some indications for primary energy and mass. 

We have used the MOCCA program (Hillas 1981) with the Sybill model (Fletcher et al. 1994) for high 
energy hadronic interactions. Twelve hundred air showers were run, at a thinning level of 10~ 6 , for two types 
of primaries (proton and iron nucleus), three primary energies (10 18 , 10 19 and 10 20 eV) and two zenithal 
angles (0° and 30°). Previous problems in the implementation of the Landau-Pomeranchuk-Migdal effect 
(Migdal 1956) and in choosing the seed of the random generator in Sybill (Prike 1998) had been corrected 
before these runs. The version of Sybill used with MOCCA was 1.5. The electromagnetic particle energy 
cutoff was set as IMeV. 
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ER and Fluctuations Results 
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In this section we show some results obtained from those twelve hundred showers, for the three quantities 
cited above. 

The ratio between the lateral densities of muon and 
electromagnetic particles is strongly dependent on 
the mass of the primary. As already shown by 
Gaisser et al (1978), the mean ratio for iron nuclei is 
60% higher than for proton primaries. At the other 
side, the applicability of this information is very de- 
pendent on the distinction between both components 
in the water Cerenkov detector. 
Here the superposition model is adopted for primary 
air nuclei interactions, and, consequently, the fluc- 
tuations in X max are also strongly dependent on the 
type of the primary particle. 

Iron-induced showers show less fluctuations in the 
depth of maximum than proton-induced ones, as can 
be seen in Figure 1. It can also be observed from 
this histogram that in about 30% of the showers 
it is impossible to distinguish between both types 

of primaries only by observing the absolute value 600 6 so 700 750 soo sso 900 

of X max . Of course the situation gets even worse 
when other primaries are included in the compari- 
son. The observed difference between mean X max 
-values for proton- and iron-induced showers, ER ■ Log(Ap e 
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Figure 1: Histograms of the X max extracted from our simulated 
data. 

g/cm 2 



is 99 ± 58 



where the error was taken 



from the widths of the gaussians fitted to the distributions. So, these fluctuations in X max wee hard to be 
beaten. 

Table 1 presents the values for ER and for the difference between the mean X max for iron and proton showers, 
observed in this sample of twelve hundred simulated showers. 
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bokolski et al. (lyyz) have concluded that the fluctuations are less sensitive to uncertainties m interaction 
models than absolute X max values themselves In their work they analyse X max and ER of events seen by 
Fly's Eye and Haverah Park experiments and claim that is impossible to distinguish between two mixtures 
of iron and proton primaries (in one case 30% protons, in the other, 80% protons), although it is possible to 
exclude pure proton or iron composition. 

3 Our Method: 

Figure 2 shows the correlation, shower by shower, between the sum of muon, electron and photon den- 
sities at 981m from shower core and its respective X max for vertical showers. Proton and iron showers 
are shown to separate themselves. Ob- 
viously, there are other nuclei induc- 
ing showers and a continuous plot is 
expected. Nevertheless, observing the 
mean behaviour of both quantities, it 
is possible to estimate, at an event to 
event basis, whether the primary is a 
light or a heavy nucleus. 
We have used the the sum of muon, 
electron (positron) and photon densi- 
ties at 981m motivated by the informa- 
tion obtained from Haverah Park ex- 
periment that one muon with energy 
higher than 400 MeV produces the 
same signal water Cerenkov tanks as 
an electromagnetic particle with en- 
ergy higher than 250 MeV, being in- 
distinguishable from it (Pryke 1996). 
We also remember that a photon pro- 
duces a pair of electron and positron 
with half of its energy each, in average. 
So, we have used these values as en- 
ergy threshold to calculate p(981m). 
The calculated p(981) will be propor- 
tional to the water Cerenkov response 
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V E M /m 2 . Similar technique was Figure 2: Plots of the correlation between p(981) and its respective X max extracted 

presented by Cortina et al (1997). The from our simulated data ' 

difference between our method and Cortina's is that they use the density of the Cerenkov light produced in the 
atmosphere and detected by the HEGRA AIROBICC experiment (Rhode et al 1996). 
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In figure 2 we can also observe that p(981) depends on the type of primary particle, p(981) for iron showers 
being higher than for proton showers at same energy. We have also studied the combined lateral distribution 

of muons, electrons and photons by fitting the function p(R) = k ■ R R o^ ■ (^-)^ to the simulated events; 
where R is the distance to the shower core in meter. This function was fitted to Haverah Park data (Pryke 
1996) for 9 < 45° and 10 17 < E < 5.10 18 eF, resulting in rj = 3.49 - 1.29 • sec{6) + 0.165 • Log{j^y), 
Ro = 4000, R± = 800 and f3 = 1.03. In the table above we present the fitted parameters to our simulated 



data. the parameter rj is equal to Z.Wd at same energy tor the Haverah Fark parametrization. 



4 Conclusions 

We have developed a method to estimate whether the primary is a light or a heavy nucleus observing 
the mean behaviour of the sum of muon, 
electron (positron) and photon densities at 
981m and its depth of the shower maxi- 
mum. The advantage of our method is that 
it allows us to estimate the type of the pri- 
mary particle at an event to event basis and 
the application of our method in the Auger 
Observatory data is very hopeful. 
We can see that the lateral distribution for 
iron showers is different from that for pro- 
ton showers. So, p(981) is also different. 
Our lateral distribution results differ from 
Hillas' results (1971) that show only a slight 
dependence on the type of the primary par- 
ticle. 

We will simulate more showers for other 
primaries and the detection of particles by 
the water Cerenkov tanks to include the ex- 
perimental efficiency. Finally, we would 
like to thank the financial support from 
FAPESP, A.M. Hillas and C.L. Pryke for 
MOCCA Code Program, and, for impor- 
tant help with technical problems, Mario 
de Castro Souza Jr.(Computing Division of Figure 3: Plots of the combined lateral distribution p versus the distance from 
IFGW - Unicamp). the shower core - 
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